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Abstract 

We show that wavefunctions in a two-dimensional (2D) electron system 
with spin-orbit coupling can be characterized by a topological quantity-the 
Chern integer due to the existence of the intrinsic Kramers degeneracy. The 
localization-delocalization transition in such a system is studied in terms 
of such a Chern number description, which reproduces the known metal- 
insulator transition point. The present work suggests a unified picture for 
various known 2D derealization phenomena based on the same topological 
characterization. 
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Metal-insulator phase transition problem in two dimensions (2D) has attracted a lot of 
attention in recent years. With time-reversal and spin rotational symmetries, all states in 
a 2D non-interacting electron system should be localized due to disorders, according to the 
Anderson localization. 1,2 In contrast to such so-called orthogonal systems, extended states 
exist in quantum Hall effect (QHE) systems 3,4 to carry current going through whole samples. 
These systems in which the time-reversal symmetry is broken by external magnetic fields are 
known as the unitary class. There are still on going debates 5-7 about whether localization- 
delocalization transition could happen in a random magnetic field case. For a time-reversal 
invariant spin-orbit coupling system, analytic and numerical studies 8-10 have been given 
consistent conclusions about the existence of a metal-insulator transition in this so-called 
symplectic class. 

These 2D quantum systems have been classified as three different universality classes, 11 
namely, orthogonal, unitary, and symplectic, in terms of symmetries of Hamiltonians. It is 
generally believed that symmetries of these systems may play important roles in determin- 
ing different localization-delocalization properties. But it is an interesting open question 
whether 2D localization-delocalization transitions in different symmetry classes are intrin- 
sically different, or, all of them (both unitary and symplectic classes) are due to the same 
universal reason? 

In the presence of magnetic fields, Thouless and co-workers 12 and others 13,14 have shown 
that the derealization property of a wavefunction is well characterized by its associated 
Chern number (which is equal to the quantized-Hall conductance in the unit of e 2 /h). As a 
topological invariant integer, Chern number represents a non-trivial topology of wavefunc- 
tion. It has been shown 14 that nodes of an eigenstate wavefunction with nonzero Chern 
number can move freely to cover the whole real space when one continuously changes the 
boundary condition. Thus a nonzero Chern number describes the extensiveness of a wave- 
function. In contrast, a state with zero Chern number should always be localized in 2D 
due to Anderson localization. 12 Numerical calculations 15 based on Chern numbers have ex- 
cellently verified localization-delocalization phase diagram in the QHE system. In random- 
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magnetic-field systems, even though the total Hall conductance becomes zero on average, 
nonzero Chern integer (or quantized Hall conductance) for single-electron states can still 
exist, and recently a metal-insulator transition has been revealed 7 numerically in terms of 
the classification of Chern integers. 

It is natural for one to ask if the derealization property in spin-orbit coupling systems 
can be similarly understood in terms of such a topological description. Without uniform or 
random magnetic field, one may expect a zero Chern number and thus a trivial topology for 
wavefunctions in their parameter space. However, this is not necessarily correct. According 
to Berry, 16 when one adiabatically changes a wavefunction in its parameter space, the Berry 
phase of the wavefunction can exhibit a monopole-like topological structure if there is a two- 
fold degeneracy. In the spin-orbit coupling system, there is a well-known two-fold Kramers 
degeneracy of eigenstates (even in the presence of random scattering) due to the time- 
reversal invariance. In this paper, we analytically show that such a degeneracy indeed 
leads to a non-trivial topological classification of wavefunctions in a spin-orbit coupling 
system which is essentially the same as in the unitary class cases, and thus one can relate 
derealization property in this system to its topological characteristics. Numerically we 
demonstrate that a topological transition happens at the point where the density of nonzero 
Chern numbers become vanishing when the system is extrapolated to infinity, and such 
a transition indeed coincides with the known critical point of metal-insulator transition 
obtained based on conventional numerical methods. 10 Therefore, a unified picture underlying 
the localization-delocalization transition in 2D systems (both unitary and symplectic classes) 
is established where all the derealizations can be essentially characterized by the same 
topological quantity- the Chern integer. 

We start with a tight-binding lattice model of noninteracting electrons with spin-orbit 
coupling. The Hamiltonian is defined as follows: 10 

H = -t Yl cfci+r -iV%2 cfa y Ci +:i: + iV ^ cfcr x c i+y + - ^ Wi&fci + H.c. (1) 

r=x,y it i i 

Here cf = (c^, cj) are fermionic creation operators, with t as the nearest-neighbor hopping 
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integral and V as hopping-spin-flipping (i.e., spin-orbit coupling) term. The strength of 
spin-orbit coupling is represented by a dimensionless parameter S — V/ \Jt 2 + V 2 . And Wi is 
a random potential with strength \wi\ < W/2 (W is adjustable). For simplicity, we assume 
no correlations among different sites for Wi (i.e., white noise limit). Such a system will be 
studied under a generalized boundary condition: + Lj) = e l9j ^/(i) with square lattice 
width Li = L 2 = L and a total number of lattice sites N = L x L (here ^ includes two spin 
components and j = 1 ,2 represent x and y directions, respectively). 

Eigenstates of H may be regarded as functions of boundary-phases 0i and 2 , and topo- 
logical property of each eigenstate (labeled by m) can be characterized in terms of relations 
with 9\ and 9 2 through the Chern integer defined as follows: 12 " 14 



in which 

F(9 1 ,9 2 ,i 



d9 1 



±J Jd9 1 d9 2 J2F(9 1 ,9 2 ;t), 



(2a) 



8i> m {9 u 9 2] i) 
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d9 2 I \ d9 2 
and the area-integral in (2a) covers a 2ix x 2tt unit cell in the #-space. The Chern number 
in (2a) can be shown 12 ' 13 to be an integer. Naively one would expect a zero Chern 
number for a system with time-reversal symmetry, since is proportional to a Hall 
conductance as mentioned earlier. In fact, one would see cancellation in (2) due to the 
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time-reversal invariance of the original Hamiltonian. However, in the spin-orbit coupling 
case, there is always a two-fold degeneracy (Kramers degeneracy) in Hamiltonian (1) at 
boundary conditions with 1; 9 2 = 7rx integer. 10 Both and its time-reversal counterpart 
state are degenerate and orthogonal here. Due to such a degeneracy, derivatives in (2b) 
would become uncertain at the degeneracy points and the contribution from these points 
needs more careful examination. 

In fact, the Chern integer is equivalent to a Berry phase (divided by 2ir) defined 
on a closed path along the boundary of a unit cell, 16 which may be converted to a surface 



integration like in (2a) by Stokes theorem. Since each degeneracy point in the parameter 
space will lead to a monopole-like singular contribution to the Berry phase, 16 a 2D #-plane 
would cut right through such a singularity point. To avoid this difficulty, one has to intro- 
duce a third parameter e to deform the surface integral around the degeneracy points. For 
example, we may add an infinitesimal term H p = — e J2i cf&zCi to the Hamiltonian. Then 
the surface integration in (2a) is understood as the one in a three-dimensional parameter 
space (#1, 9 2 , e) (here e <C max.(t,V,Wi)/N so that it has no effect when one is sufficiently 
far away from the degeneracy points). Due to the cancellation in the regions far away from 
the degeneracy points, Eq.(2a) is then reduced to 



in which Sk (k — 1 — 4) represent four small surface-integrals covering the degeneracy points 
(0,0), (ir, 0), (0, 7r), and (jr, n) within a 27r x 27r unit cell. The minimum size of these areas 
is decided by the strength e and the final result should not depend on e at e — > + . One may 
find that this regularization procedure resembles the one in a magnetic phase transition 
where an infinitesimal external symmetry-breaking field is usually introduced to lift the 
ground-state degeneracy for a spontaneous magnetization. Here H p plays a role to lift the 
Kramers degeneracy such that the hidden topological characterization becomes explicitly 
shown (see below). 

Denote each degeneracy point as (0° fe , 9^ k ) (k — 1 — 4) and suppose that eigenstates 
are already known at these points for a given e. Then one can determine eigenstates at an 
arbitrary boundary phase (9±, 9 2 ), which is close enough to (0° fe , 9®k), by a perturbative 
way. The following procedure is similar to that in Ref. 17 in treating the level-crossing 
problem. The change of the boundary condition may be incorporated into the Hamiltonian 
as a perturbation, and in terms of the complete set of basis obtained at (9^,9^), such a 
change in Hamiltonian can be generally written as 
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where 2x2 matrix elements are related to the (2n-l)-th and (2n)-th eigenstates (assuming 
an energy-increase order for all the eigenstates), which correspond to a pair of two-fold 
degenerate eigenstates of H. The off-diagonal terms may be expanded to a leading order of 
(9i — 9\ k ) and (9 2 — 9 2k ) (If the leading term happens to be zero, higher order expansions 
in AH(9i,9 2 ) will be needed, and integer Chern numbers can be similarly obtained 17 ). The 
general form is A n (9 ± , 9 2 ) = - 9° lk ) + (3 k (9 2 - 9° 2k ) with a k = ^ < * 2 n-i(^ fe , Q 2k )\H + 
H p \* 2n (9° lk ,9° 2k ) > and fa = ^ < ^ 2n -i(9° lk , 9° 2k )\H + H p \V 2n (9° lk , 9° 2k ) >. By diagonalizing 
the 2x2 matrix of the Hamiltonian for each n, we get the (2n-l)-th and (2n)-th eigenstates at 
the boundary phase (9±, 9 2 ) sufficiently close to (9® k , 9 2k ). Then according to (3), the Chern 
numbers for the (2n-l)-th and (2n)-th eigenstates are given by the following integers: 17 

^^EsgnNa^-)]; (4a) 

C (2n) = _ C {2n-l)^ ^ 

whose values will not depend on e = + . (Here if Im(a*/3j) = 0, we define sgn[lm(a*/3j)} = 
in (4a).) Therefore, eigenstates in the spin-orbit coupling system can be indeed classified 
in terms of nontrivial Chern integers, and two degenerate eigenstates always have Chern 
numbers with opposite signs such that the total sum of them (and the sum of the Hall 
conductances) is still zero in consistence with the time-reversal symmetry. Here each de- 
generacy point gives rise to a contribution similar to a magnetic monopole in unit strength 
1/2 as originally discussed by Berry. 16 Nonzero Chern integers have been used to character- 
ized extended states in both the QHE and random-magnetic- field systems. 14 ' 15 ' 7 ' 18 By the 
same token, the derealization problem in the spin-orbit coupling system should be simi- 
larly described here, after the Chern-number classification of wavefunctions is established. 
In the following, we shall use this topological characteristics to determine the localization- 
delocalization transition in the present system. 

We define a density of states p ex t(£, N) for the eigenstates with nonzero Chern number at 
energy e. p ex t is related to the density of states for extended states. 15 ' 7 Exact diagonalization 
of H is carried out at four boundary conditions where the two-fold degeneracy happens, and 
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then Chern numbers are determined through the formula (4). The total density of states 
p{e, N) and the extended one p ext (s, N) are calculated as a function of lattice size N (iV=36, 
64, 144, 256 and 576), which is averaged over random potential configurations (100 — 5000 
random configurations depending on sample sizes) and also averaged over a small energy- 
width As = 0Ay/V 2 + t 2 in the neighborhood of energy e to reduce statistics error. The 
ratio p ext (e,N)/p(e,N) at the band center e = is presented in Fig. 1 with a log-log 
plot. Different curves correspond to different random strength Ws. The spin-orbit coupling 
strength S has been chosen to be S = 0.5 with \/V 2 + t 2 = 1 (in order to compare with 
the results in Ref. 10). The ratio p ex t{s = 0,N)/p(e = 0,N) monotonically decreases 
as parameter W increases, which reflects the fact that eigenstates become more and more 
localized at stronger disorders. For smaller W (W=3, 4 and 5 as shown in the Fig. 1), 
the ratio p ex t{s = 0,N)/p(e = 0,N) shows a flat curve with a slight increase at larger 
sample sizes, suggesting that a finite value of p ext {e = 0,N)/p(e = 0,N) may be reached 
if it is extrapolated to infinity lattice size. The trend becomes opposite at larger strengths 
of random potential, where the ratio p ex t(£ = 0,N)/p(e = 0,N) becomes monotonically 
decreasing with larger sample sizes. In fact , all the data at W > 5 shown in Fig. 1 can be 
nicely fit into a straight line, which means p ext (e = 0, N)/p(e = 0, N) oc N~ x (with x > to 
be given below). These finite-size scalings show that at a sufficiently large W, the density of 
states with nonzero Chern number is always extrapolated to zero at large lattices. Similar to 
the cases in the QHE and random-magnetic-field systems, it means that all the eigenstates 
here have trivial topology and should be localized. 

The exponent function of the disorder strength W is shown in Fig. 2. It is 

interesting to see that all data fall on a straight line. Thus one can determine the critical 
strength of disorder W c at x c = by extrapolation in Fig. 2. Here W c decides a topological 
phase transition point. At W > W C) p ext ~ N~ x — > at N — > oo, while at W < W C) 
Pext — ► a finite constant at large lattice size. As have been known in the QHE and random- 
magnetic-field systems, such a topological transition physically corresponds to a localization- 
delocalization transition. In fact, W c determined from Fig. 2 is approximately 5.62 ± 0.08, 
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which is quite close to the value W c = 5.74 obtained by conventional numerical methods 10 
with much larger lattice sizes ~ 64 x 10 5 (A slightly smaller W c in our result may be due 
to the fact that a finite energy average around e = is made in calculation which could 
effectively lower W c ). 

The above analytic analysis and numerical calculation have shown that derealization 
occurs in the spin-orbit coupling case is basically the same as that in the unitary cases. It 
suggests that all localization-delocalization transitions in 2D non-interacting systems may 
be universally described as a topological phase transition in terms of the Chern integer. 
The conventional understanding of metal-insulator transition based on symmetry classes of 
Hamiltonians may not be essential. For instance, by applying a small but finite magnetic 
field one may change a spin-orbit coupling system from symplectic class to unitary class by 
destroying the time-reversal symmetry. We have checked this case numerically and found 19 
no sudden change in the distribution of Chern integers (thus the topological transition), 
which actually behaves quite smoothly with small magnetic fields. This fact has already 
been noticed before in the behavior of localization length. 20 The universality has also been 
shown in the study of the critical behavior of metal-insulator transition, where it has been 
found 21 that correlation dimensions of spectral measure D 2 and fractal eigenstate D 2 as well 
as exponent rj describing the energy correlations of the critical eigenstates in spin-orbit cou- 
pling system are very close to those obtained for the QHE systems. 22 Experimentally, both in 
random magnetic field 23 and spin-orbit coupling 24,25 systems, magnetoconductances are neg- 
ative in the delocalized region ( here "derealization" means the localization length is at least 
larger than the sample size), while positive in the localized region. This common feature can 
be understood as that in the delocalized region, a finite fraction of eigenstates carry nonzero 
Chern number and magnetic fields suppress original phase coherence and reduce the number 
of states having nonzero Chern numbers. 19 Thus the magnetic field plays a role in reducing 
derealization in the "delocalized" region and leads to negative magnetoconductance. In the 
opposite case, where original wavefunctions are topology trivial, magnetic field will be in 
favor of derealization to give rise to a positive magnetoconductance. 26 This explains why 
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metal-insulator transition is generally accompanied by a sign change of megnetoconductance 
experimentally in spite of different symmetry classes. 

In conclusion, the localization-delocalization transition in a spin-orbit coupling system is 
studied from a topological point of view. It is shown that wavefunctions in this time-reversal 
invariant system can be characterized by Chern integers as a result of the intrinsic Kramers 
degeneracy. We have determined the critical strength of disorder W c for the metal-insulator 
transition point at energy e = based on such a Chern-integer characterization, which is in 
excellent agreement with earlier numerical works 10 by localization-length scalings. Therefore, 
various known localization-delocalization transitions in 2D non-interacting systems now can 
be all described in terms of a topological characterization. Whether there is a far-reaching 
physical implication behind it needs to be further studied. 
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FIGURE CAPTION 



Fig. 1. The ratio of the density of extended states over the total density of states is 
plotted as a function of sample size N at energy e — 0. W is the strength of random scat- 
tering. Error bars of the data for W < 5.0 are in the same order of magnitude as those for 
W > 5.0 but are not marked in the figure for clarity. 

Fig. 2. The exponent x vs. the strength of random scattering W. The data is extrapo- 
lated to x — 0, giving a critical strength W c = 5.62 ± 0.08. 
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